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PREFACE 

The  work  reported  herein  was  conducted  by  the  Arnold  Engineering  Development 
Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at  the  request  of  the  Air  Force 
Engineering  and  Service  Center,  Tyndall  Air  Force  Base,  Florida.  The  results  of  the  research 
were  obtained  by  ARO,  Inc.,  AEDC  Division  (a  Sverdrup  Corporation  Company), 
operating  contractor  for  the  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee,  under 
ARO  Project  Numbers  E34P-R4A  and  P34I-02A.  The  data  analysis  was  completed  on 
June  4,  1979,  and  the  manuscript  was  submitted  for  publication  on  August  20,  1979.  The  Air 
Force  project  managers  were  Capt.  Harvey  Clewell,  of  Tyndall  Air  Force  Base,  and  E.  L. 
Hively,  Directorate  of  Technology,  AEDC.  The  experimental  phase  was  conducted  at  the 
United  Technologies  Research  Center  Facility  in  East  Hartford,  Connecticut. 

The  authors  wish  to  express  their  appreciation  to  M.  F.  Zabielski,  L.  G.  Dodge,  and  M. 
Colket,  UTRC,  for  their  cooperation  and  dedication  during  the  experimental  phase  of  this 
program.  In  addition,  the  authors  gratefully  acknowledge  the  assistance  provided  by  W.  E. 
Myers,  ARO,  Inc.,  from  system  buildup  through  data  acquisition. 
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1.0  INTRODUCTION 

Differences  between  measurements  of  nitric  oxide  (NO)  emissions  from  combustion 
sources  using  conventional  sample  extractor  techniques  and  those  using  an  in  situ  ultraviolet 
(UV)-absorption  (optical)  technique  have  been  reported  (Refs.  1,  2,  and  3).  The  Federal 
Aviation  Administration  (FAA)  is  directing  an  Interagency-sponsored  program  to  determine 
the  cause  of  the  differences  found  between  the  measurements  using  the  two  techniques.  As 
part  of  this  investigation,  United  Technology  Research  Center  (UTRC)  has  provided  several 
combustion  sources  to  evaluate  the  differences  between  the  optical  and  probe  techniques 
(Ref.  4).  Arnold  Engineering  Development  Center  (AEDC)  personnel  have  made 
measurements  using  the  UV-absorption  technique  (developed  at  AEDC)  at  the  UTRC 
facility  in  East  Flartford,  Connecticut,  and  the  results  of  the  first  phase  of  that  study  are 
reported. 

The  Interagency  NO  Investigation  has  three  major  tasks.  Task  1  is  the  development  of 
calibration  techniques  for  the  optical  measurement  systems;  Task  2  emphasizes  determining 
the  relevant  chemistry  occurring  in  the  sample  extraction  and  in  the  transfer  lines  to  the  NO 
measurement  instrumentation;  Task  3  includes  both  optical  and  probe  sampling  on  three 
progressively  more  complicated  combustion  systems.  UV-absorption  measurements  will  also 
be  made  by  UTRC  using  basically  the  same  approach  as  that  used  by  AEDC,  and  an 
infrared  correlation  technique  is  also  used  (Ref.  4).  The  results  of  all  measurements,  both 
probe  and  optical,  will  be  compared  and  analyzed. 

This  report  summarizes  the  results  of  measurements  made  by  AEDC  during  the  optical 
calibration  phase  (Task  1)  of  the  overall  program.  During  the  Task  1  phase,  some 
discrepancies  were  identified  by  UTRC  personnel  (Ref.  4)  in  the  original  theoretical  optical- 
transmission  model  developed  at  AEDC  for  nitric  oxide  concentration  measurements.  The 
model  has  since  been  corrected  by  AEDC/ARO  personnel;  a  review  of  the  theory  is  included 
in  this  report.  (UTRC  has  also  developed  its  version  of  the  model,  which  is  described  in  Ref. 
4.)  Model  predictions  are  compared  to  data  points  obtained  from  the  UTRC  calibration  cell, 
flowing-gas  heater,  and  flat-flame  burner.  These  data  points  span  the  temperature  range 
from  300  to  1,750  K  at  a  pressure  of  1  atm.  Agreement  of  the  optical  measurements  with 
calculated  inputs  of  NO  concentration  in  conjunction  with  probe  profiles  is  generally  within 
the  uncertainty  of  the  probe  measurements. 
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2.0  REVIEW  OF  MODEL 


2.1  THEORY 

A  theoretical  model  for  the  transmission  of  Doppler-broadened  resonance  line  radiation 
through  absorbing  media  with  combined  Doppler  and  pressure  broadening  was  given  by 
Davis,  McGregor,  Few,  and  Glassman  (Ref.  5)  with  application  to  the  y  band  of  NO.  As 
shown  by  Ref.  4,  the  formulation  given  in  Ref.  5  was  incorrect  in  that  Hund’s  case  (b) 
statistics  (Ref.  6)  were  used,  the  spin  splitting  was  not  properly  taken  into  account,  and  the 
normalization  of  the  line  oscillator  strength  was  inconsistent  with  the  literature  (Refs.  7  and 
8).  The  method  of  calculation  of  the  partition  function  in  the  model  has  now  been  changed 
to  conform  to  the  angular  momentum  coupling  intermediate  between  Hund’s  cases  (a)  and 
(b)  in  accordance  with  Keefer  (Ref.  9)  and  as  discussed  in  Ref.  4.  A  brief  review  of  the  theory 
is  presented  in  the  following. 

The  transmission,  f  ,  of  an  isolated,  Doppler-broadened  emission  line  through  a 
homogeneous  medium  of  length,  f,  is  given  by 


where  I^P is  the  intensity  of  the  source  line  at  line  center,  r-1  is  the  frequency  at  line  center, 
(Asrj)D  is  the  half-width  of  the  Doppler-broadened  line,  and  k(j>)  is  the  monochromatic 
absorption  coefficient  of  the  medium.  The  half-width  of  the  source  line  depends  on  the 
source  temperature  and  is  given  by 


(2) 


where  Ms  is  the  molecular  weight  of  the  emitting  molecule,  Ts  is  the  source  temperature,  x  is 
Boltzmann’s  constant,  and  c  is  the  speed  of  light. 
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The  absorption  coefficient,  k(^),  contains  contributions,  k,(^),  from  all  lines  near  the 
frequency,  //’■  Each  of  these  lines  will  be  broadened  both  by  the  Doppler  effect  and  by 
collisions.  The  absorption  coefficient  of  the  ith  line  is  given  by  the  Voigt  profile, 


fie-y 


a; 


("i-T 


dy 


(3) 


where  kf  is  the  centerline  Doppler  absorption  coefficient  of  the  ith  line,  and  the  broadening 
parameter,  Cf,  is  a  parameter  proportional  to  the  ratio  of  the  collision-broadened  half-width, 
(^a^i)L,  to  the  Doppler-broadened  half-width,  (Aai/j)D,  in  the  absorbing  medium.  The 
broadening  parameter,  0,  equation  is  as  follows: 


(4) 


The  parameter  cof  in  Eq.  (3)  is  the  Doppler  function 


2{v-vty 

(^a^i)  D 


y/  In  2 


(5) 


It  can  be  shown  (Appendix  A  and  Ref.  10)  that  for  binary  collisions 


a  =  cp/T 


(6) 


where  p  is  the  pressure,  T  is  the  temperature,  and  C  is  a  constant  which,  for  any  given 
mixture  of  broadening  species,  b,  of  mass,  Mb,  mole  fraction,  fb,  and  broadening-collision 
cross  section,  crab,  is  given  by 
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c  = 


A 


3/2 


£  fbaab 
b 


(7) 


where  X  is  the  wavelength  and  dab  is  the  optical  diameter.  The  appropriate  absorption 
coefficient  in  Eq.  (1)  must  include  the  contributions  due  to  all  absorbing  lines  that  overlap 
the  source  line,  j;  thus, 


k(v)  = 


Z  k,M  -Erf 

[  J  S#— CK 


fi2 


dy 


(8) 


There  may  be  several  emission  lines  within  a  given  instrument  bandpass,  Av{.  The 
contributions  of  each  of  these  lines  must  be  summed  to  obtain  the  transmission  within  the 
given  instrument  bandpass.  Thus,  the  transmissivity  within  the  bandpass,  Aiq,  is  given  by 
substitution  of  Eq.  (8)  into  Eq.  (1),  summing  over  the  emission  lines,  j,  and  dividing  by  the 
total  intensity  of  the  source  within  the  bandpass  as  follows: 


where  k(j')  is  given  by  Eq.  (8). 

The  line  center  absorption  coefficient  for  a  single  transition  out  of  the  12  possible 
transition  from  a  given  J"  is  given  (Ref.  5)  by 
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where  i  refers  to  the  particular  n'v' J'K'n"v"J"K"  transition  in  the  A2£  — x2n  7-band 
system  of  NO.  The  constants  in  parenthesis  are  the  electron  charge,  e,  the  electron  mass,  Me, 
and  the  speed  of  light,  c,.  The  population  of  the  v"  vibrational  level  in  the  J"K"  rotational 
sublevel  and  one  of  the  parity  levels,  p,  ( +  or  -)  of  the  2Hn  state  in  the  intermediate  angular 
momentum  coupling  case  between  Hund’s  cases  (a)  and  (b)  is  given  (Ref.  9)  by 


N(v  "  J"K"mp) 
^total 


j 

-  G(v" )  +  F  (J") 

m 

I 

(2 J"  +  1 )-  exp  / 

kT 

S 

2  Y. 

|  exp 

2 

r-G(v)i  £ 

00 

£  (2 J  +  1)  exp 

J=0 

m. 

V 

L  kT  J  m=l 

^  kT  j( 

(ID 


where  T  is  the  temperature,  G(v")  is  the  vibrational  energy  (cnr1).  Fm(J")  is  the  rotational 
energy  (cm1)  for  either  the  2H]/2  state  (m  =  1)  or  the  2n3/2  state  (m  =  2)  and  for  purposes 
of  the  population  calculation  is  approximated  by  the  Hill  and  Van  Vleck  formulas  (Ref.  6) 


(J  +  V2)2 

J4(j  +  %f  +  y(y  -  4; 

>  -D,J4 

(12a) 

(j  +  V2)  2  -  1  +  1/2 

%)2  +  y(y  -  4) 

-dT  +  i)4 

(12b) 

where  Y  =  A/Bv  where  A  is  the  splitting  constant,  Bv  is  the  inertial  rotational  constant,  and 
Dv  is  the  contrifugal  rotational  constant.  The  oscillator  strength  for  the  transition  is  given  by 


f(n’v’  J’K*  h"v" J"K") 


sG*  K1  J" K")  jyG'K*  J"K")  V,  A 

2J"  +  1  v(v'  v"  ) 


(13) 


where  p(J'K' J"K")  and  Kv'v'')  are  the  frequencies  of  the  line  and  band  centers  and  where 
the  normalization  for  the  Honl-London  factor  S(J'K'J"K")  and  the  normalization  of  the 
band  oscillator  strength  follow  that  of  Tatum  (Ref.  7)  and  are  consistent  with  Pery-Thorne 
and  Banfield  (Ref.  8).  The  Honl-London  factors  are  those  of  Earls  that  had  been  defined 
earlier  (Ref.  5). 


9 


A  EDC-TR-79-65 


The  line  positions,  v®,  previously  used  were  those  of  Deezi  (Ref.  11).  Improved  data  have 
recently  become  available  (Ref.  12),  and  these  data  have  been  correlated  by  UTRC  (Ref.  4) 
into  a  consistent  formulation.  These  new  line  positions  have  been  incorporated  into  the 
AEDC  model. 

The  source  parameters  are  the  peak  line  intensity,  I°9,  and  the  Doppler  half-width  of  the 
lines,  (Asj'j)D.  It  is  sufficient  to  determine  only  the  relative  peak  line  intensity  since  any 
absolute  magnitude  will  cancel  in  Eq.  9.  The  relative  peak  line  intensity  can  be  determined 
by  direct  measurement  for  resolved  lined;  the  intensities  of  the  unresolved  lines  can  be 
determined  from  the  upper  state  population  implied  by  the  measured  line  intensities. 
Elowever,  for  low-pressure,  Doppler-broadened  source  lines  in  the  ultraviolet,  direct 
measurement  of  either  the  half-width  or  the  source  temperature  is  difficult.  Section  2.2 
describes  an  indirect  method  of  determining  source  temperature. 

The  absorber  parameters  are  the  absorption  coefficient  at  line  center,  kf,  and  the 
quantity  d,  related  to  the  collisional  broadening.  If  the  temperature  and  density  within  the 
absorption  path  are  known,  then  k° is  determined  by  Eq.  (1 1)  if  the  band  oscillator  strength, 
f(v'v"),  is  known.  Thus,  when  kf  is  known,  the  broadening  coefficient,  C,  of  Eq.  (6)  can  be 
determined  from  the  value  of  broadening  parameter,  ff,  required  to  obtain  a  match  between 
the  measured  transmissivity  and  the  transmissivity  calculated  from  the  model  based  on  Eqs. 
(9)  through  (13). 

Successful  application  of  the  model  represented  by  Eqs.  (9)  through  (13)  requires  a 
known  source  line  intensity  distribution,  source  temperature,  band  oscillator  strength,  and 
broadening  coefficient. 

2.2  PARAMETERS 

The  model  represented  by  Eqs.  (9)  through  (13),  together  with  carefully  controlled 
laboratory  experiments,  can  be  used  to  determine  the  source  temperature,  Ts,  and  the 
broadening  parameter,  (f,  when  the  source  line  positions  and  intensities  are  specified  and  the 
band  oscillator  strength  is  known.  Careful  measurements  of  oscillator  strength  have  been 
made  for  the  NO  molecule  using  several  methods  (e.g.,  Refs.  8,  13,  and  14).  Review  of  these 
measurements  shows  the  value  of  4.09  ±  0.1  x  104  of  Hassan  et  al.  (Ref.  13)  for  the 
oscillator  strength  of  NO  to  be  the  most  reliable. 

2.2.1  Source  Lines 

The  line  positions  of  the  (0,0)  band  of  the  y  system  of  NO  as  specified  by  Ref.  4, 
interpreting  Engleman’s  data  (Ref.  12),  have  been  adopted  as  the  most  accurate  available. 
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The  source  line  intensity  distribution  of  the  AEDC  capillary  discharge  lamp  has  been 
remeasured  with  resolution  sufficient  to  identify  many  lines  lying  in  each  of  the  independent 
branches  of  the  (0,0)  band.  No  changes  from  the  previously  reported  intensity  distribution 
(Ref.  5)  were  found.  It  should  be  noted  that  the  line  discharge  lamp  and  its  operating 
conditions  (pressure,  gas  mixture,  and  voltage  applied)  and  different  designs  (e.g.,  hollow 
cathode  configuration  as  used  in  Ref.  4)  would  be  expected  to  exhibit  a  different 
distribution.  Furthermore,  there  is  no  reason  to  expect  the  intensity  distribution  to  follow  a 
Boltzmann  distribution  of  the  population  because  of  the  complex  excitation  mechanism 
within  the  electrical  discharge. 

2.2.2  Source  Temperature 

Determination  of  the  effective  source  temperature  hinges  on  the  fact  that  at  low  pressure 
the  transmissivity  becomes  independent  of  the  other  unknown  parameter,  the  collision¬ 
broadening  parameter.  In  the  low-pressure  limit,  the  combined  Doppler-  and  collision- 
broadened  absorption  profile  [Eq.  (3)]  reduces  to  a  Gaussian  profile  whose  half-width 
depends  only  upon  the  absorber  temperature.  Thus,  if  the  absorber  temperature  and  number 
density  are  known,  and  the  measured  value  of  oscillator  strength  is  used,  the  only  unknown 
in  Eq.  (9)  is  the  Doppler  half-width  of  the  source.  This  half-width  is  directly  related  to  the 
source  temperature  through  Eq.  (2). 

The  experimental  apparatus  described  in  Ref.  5  was  used  to  make  a  carefully  controlled 
set  of  transmission  measurements.  A  sequence  of  transmissivity  measurements  at  the  second 
bandhead  was  made  over  a  range  of  pressures,  with  the  NO  number  density  and  temperature 
held  constant.  In  order  to  determine  the  source  temperature,  the  transmissivity  for  a  given 
number  density  and  temperature  was  determined  from  the  intercept  of  the  transmissivity 
versus  pressure  curve  (Fig.  1)  in  the  limit  of  zero  pressure  (corresponding  to  (3  =  0).  This 
value  of  transmissivity  was  used  together  with  Eq.  (9)  in  the  limit  where  (3  approaches  zero  to 
determine  the  source  temperature  by  adjusting  source  temperature  to  force  a  match  between 
the  calculated  and  extrapolated,  measured  transmissivities.  A  value  of  source  temperature  of 
950  ±  25  K  was  obtained  by  this  exercise.  Again,  it  should  be  noted  that  the  source 
temperature  is  a  characteristic  of  the  lamp  design  and  is  expected  to  vary  for  different 
designs. 

2.2.3  Broadening  Coefficient 

Once  the  source  temperature  has  been  determined,  the  only  unresolved  parameter  in  Eq. 
(9)  is  (3,  the  parameter  proportional  to  the  ratio  of  collision-  to  Doppler-broadening  half¬ 
widths. 
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p/T,  atm/K 

Figure  1,  Least-squares  fit  to  measured  transmissivity  versus  p/T. 

A  series  of  experiments  was  performed  for  known  concentrations  of  NO  and  various 
diluent  gases  over  a  range  of  pressures  from  15  to  760  torr.  The  measured  values  of 
transmissivity — together  with  Eqs.  (9)  through  (13) — were  used  to  determine  the  value  of  a 
for  each  experimental  condition.  A  least-squares  curve  fit  to  these  data  was  used  to 
determine  the  value  of  C  in  Eq.  (6),  and  the  broadening  cross  section  was  determined  from 
Eq.  (7).  Since  the  calibration  gas  was  a  mixture  of  NO  in  N2)  the  coefficient,  C,  for  both  Ar 
and  Co2  was  determined  by  use  of  Eq.  (7).  Figure  2  is  a  plot  of  Q  versus  p/T  for  some  100 
data  points  using  N2  as  the  broadening  gas  and  showing  the  least-squares  straight  line  fit  to 
the  data.  The  results  for  the  broadening  of  No  by  N2,  A,  and  C02  are  given  as  follows: 

Collision  0,-7-  Cross  Section,  Collision  Diameter, 

Partner  atm  IQ-16  cm2  A° 

N2  2,390  ±  220  556  ±  53  13,3  ±  0,6 

A  2,150  ±  200  547  ±  53  13.2  ±  0.6 

C02  2,090  ±  200  540  ±  53  13.1  ±  0.6 
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It  should  be  noted  that  at  the  higher  pressures  where  the  absorbing  lines  are  considerably 
broader  than  the  source  lines,  the  transmissivity  becomes  relatively  independent  of  the 
source  temperature  determined  from  the  low-pressure  measurements.  However,  at 
intermediate  pressures,  where  source  line  width  and  absorber  line  width  are  comparable,  the 
correct  source  temperature  is  essential  for  an  accurate  determination  of  the  collision¬ 
broadening  parameter  by  the  method  used  here. 

The  magnitude  and  the  near-constant  value  of  the  collisional  cross  section  for 
broadening  of  NO  7-band  lines  by  the  three  molecules  studied  suggest  that  long-range 

electronic  forces  are  responsible  for  the  broadening  and  that  motion  of  the  nuclei  is  not 
affected.  For  data  reduction  in  the  UTRC  experiments,  the  following  equation  was  used 


p/I,  atm/K 

Figure  2.  Values  of  the  spectral  broadening  parameter,  S  ,  as  a 
function  of  p/T  for  the  (0,0)  7  band  of  NO. 
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where  C0  =  1,656  K/atm  from  the  measurements  reported  here.  It  should  be  emphasized 
that  Eq.  (14)  when  applied  to  other  gases  represents  an  assumption.  For  the  UTRC 
experiments,  gases  of  interest  are  N2,  A,  H2o,  and  02.  Thus  C  for  H20  and  02  are  assumed 
to  follow  Eq.  (14)  purposes  of  data  reduction. 

2.2.4  Model  Reproducibility 


If  the  source  line  intensity  distribution,  source  temperature,  oscillator  strength,  and 
collision  diameter  have  all  been  correctly  determined,  the  model  should  accurately  predict 
the  transmissivity  over  the  complete  pressure  and  temperature  range  covered  by  the 
measurements.  In  Fig.  3  the  model  represented  by  Eqs.  (9)  through  (13),  using  the 
parameters  determined  as  described  above,  is  compared  to  the  laboratory  transmissivity 
measurements  at  the  second  bandhead  over  the  pressure  range  from  15  to  760  torr  for  an 
absorption  cell  density  of  2. 14  x  1014  cm-3  at  room  temperature.  The  model  was  also  checked 
successfully  using  a  wide  range  of  other  absorption  cell  data  taken  before  and  after  the 
UTRC  measurement  program.  The  standard  deviation  for  all  points  of  Fig.  3  is  ±0.5 
percent.  This  agreement  between  predicted  and  measured  transmissivity  over  a  large 
pressure  range  indicates  the  accuracy  of  the  current  AEDC  model  and  the  correct 
determination  of  the  model  parameters  within  the  pressure  and  temperature  ranges  used  in 
the  laboratory.  It  is  also  noteworthy  that  the  model  now  predicts  the  first  bandhead 
transmissivity  accurately,  contrary  to  the  previous  model. 
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Figure  3.  Experimental  measured  transmissivity  compared  to  model 
predicted  transmissivity. 
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3.0  OPTICAL  CALIBRATION  RESULTS 

3.1  UTRC  CALIBRATION  SOURCES 

The  UTRC  high-temperature  calibration  system  (flowing-gas  heater  and  flat-flame 
burner)  is  described  fully  in  Ref.  4.  It  was  designed  to  minimize  both  gas  phase  and  surface 
reactions  by  using  argon  as  the  carrier  gas  in  the  heater  system  and  as  a  diluent  in  the  lean 
H2/02  flat-flame  burner.  Also,  under  these  conditions  the  loss  of  NO  in  probes  should  be 
minimized  because  of  the  lack  of  hydrocarbons  and  free  hydrogen  in  the  high-temperature 
gas  flow,  the  agents  suspected  to  be  responsible  for  the  NO-destroying  probe  reactions  (Ref. 
15).  Therefore,  comparison  between  probe  and  optical  results  should  provide  a  reasonable 
calibration  for  the  optical  technique.  A  static  absorption  cell  (Ref.  4)  was  also  used  to 
provide  a  datum  for  optical  measurements  by  all  parties.  These  three  sources — static 
absorption  cell,  flowing-gas  heater  and  flat-flame  burner,  and  optical  system — will  be 
described  briefly  in  this  report. 

3.1.1  Static  Absorption  Cell 

The  static  cell  was  a  stainless  steel  pipe  18.6  cm  in  length  by  2.03  cm  in  diameter  with 
quartz  windows  at  either  end.  A  gas  mixing  and  filling  system  described  in  Ref.  4  permitted 
evacuation  of  the  cell  and  filling  with  calibration  gas  to  a  known  pressure  and  NO 
concentration.  The  cell  was  operated  at  room  temperature. 

3.1.2  High-Temperature  Calibration  Facility 

The  high-temperature  calibration  facility  housed  either  the  flowing-gas  heater  (FGH)  for 
temperatures  up  to  1,000  K,  or  the  flat-flame  burner  (FFB)  for  temperatures  up  to  2,000  K. 
A  rectangular,  water-cooled,  stainless-steel  shroud  22  cm  by  14  cm  and  26  cm  high  prevented 
atmospheric  entrainment  and  contained  the  hot  gas  flow  which  was  exhausted  through  a 
vent  system  to  the  outside  of  the  building.  The  FGH  or  the  FFB  sources  were  mounted  inside 
the  shroud  at  the  base  of  the  rectangular  duct.  A  gas-mixing  and  metering  system  supplied 
the  desired  gases  to  the  sources.  Viewing  ports,  3,4  cm  in  internal  diameter  (ID)  by  10.4  cm 
in  length,  with  optical  windows  mounted  at  either  end  provided  a  viewing  path  through  the 
center  of  the  shroud  along  the  22-crn  path.  After  the  problem  of  calibration  gas  buildup 
within  the  view  tubes  was  encountered,  a  purge  of  dry  nitrogen  was  introduced  just  inside 
the  optical  windows. 

Flowing-Gas  Heater — A  quartz  bead  bed  is  heated  through  radiation  from  electrically 
powered  ceramic  heating  units  to  the  desired  temperature.  Metered  calibration  gas  diluted 
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with  either  N2  or  A  flows  through  the  quartz  bed  and  is  heated  by  contact  with  the  quartz 
surfaces.  The  calibration  gas  flows  through  a  17.4-  by  9.2-cm  (158.3-cm2)-area  quartz  duct 
into  the  rectangular  shroud.  The  calibration  gas  is  surrounded  by  a  narrow  sheath  of  buffer 
gas,  N2,  which  has  a  temperature  near  that  of  the  calibration  gas.  The  stream  overall  NO 
concentration  is  controlled  by  the  metered  flow.  The  concentration  and  temperature. profiles 
are  influenced  by  the  optical  port  purge  gas,  by  normal  mixing  with  the  external  stream,  and 
by  the  heater  characteristics  (Ref.  4).  The  stream  temperature  and  concentration  profiles 
along  the  optical  path  were  measured  by  UTRC  using  sampling  and  thermocouple  probes 
located  a  distance  of  4.0  cm  above  the  heater  exit  plane.  The  optical  axis  was  also  located  at 
a  distance  of  4.0  cm  above  the  exit  plane. 

Sampling  probes  of  various  designs  were  used  to  access  different  portions  of  the  optical 
path.  (The  standard  probe  was  too  large  for  the  area  near  the  stream  edges.)  The  probes  are 
described  in  Ref.  4.  Sample  line  measurements  were  made  using  either  a  chemiluminescent 
analyzer  or  a  mass  spectrometer,  the  majority  of  the  measurements  made  by  the  latter. 
Concentration  profiles  provided  by  UTRC  to  correspond  to  the  optical  measurements  are 
shown  in  Figs.  4  and  5.  The  profiles  are  identified  by  profile  designation  in  Table  1  with  a 
particular  optical  measurement.  The  profiles  are  shown  normalized  to  values  determined  by 
the  gas-metering  system.  The  available  probe  measurements  for  the  centerline  are  also  given 
in  Table  1. 


Probe  Profile 
Sym  Design  Designation 


Figure  4.  Normalized  NO  concentration  profiles  of  the  flowing-gas 
heater  at  room  temperature. 
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Distance  from  Centerline,  cm 


Figure  5.  Normalized  NO  concentration  profiles  of  the  flowing-gas 
heater  at  elevated  temperature. 

The  gas  temperature  profiles  along  the  optical  path  for  the  FGH  were  obtained  by  a 
Chromel®-Alumel®  thermocouple  probe  (O.OlO-in.-diam  wire  and  0. 015-in. -diam  junction). 
These  profiles  are  shown  on  Fig.  6  for  the  data  points  indicated  in  Table  1.  Radiation 
corrections  were  not  made  to  these  data  because  they  were  thought  to  be  negligible  (Ref.  4). 
The  nominal  temperature  for  each  data  point  is  given  in  Table  1.  The  nonuniformity  in  the 
temperature  profile  has  been  attributed  by  UTRC  to  transient  conditions  in  the  heater;  these 
conditions  are  peculiar  to  the  particular  data  points  listed. 


Table  1.  Operating  Conditions  for  Flowing-Gas  Heater 


Data 

Temperature 

Average 

Concentration 

Centerline 
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Temperature, 
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Distance  from  Centerline,  cm 

Figure  6.  Temperature  profiles  of  the  flowing-gas  heater 

Flat-Flame  Burner — The  flat-flame  burner  consisted  of  a  lean  burning  H2/02  premixed 
porous-plate  burner  in  which  NO  and  A  calibration  gases  were  mixed.  Mixture  contents  in 
the  form  of  mole  fractions  are  given  for  each  optical  measurement  point  in  Table  2.  The 
interior  of  the  rectangular  burner  measures  17.5  by  9.2  cm  (161  cm2).  A  surrounding  buffer 
gas  (area  76  cm2)  of  lean  burning  FF/CF  and  argon  was  also  supplied. 

The  gas  sampling  of  the  flat-flame  burner  at  a  distance  of  from  1  to  2  cm  (depending  on 
the  data  point)  from  the  burner  face  was  accomplished  in  the  main  part  of  the  flow  by  water- 
cooled  quartz  probes.  Near  the  wall  (less  than  1.5  cm)  an  uncooled  quartz  probe  was  used. 
Concentration  profiles  for  the  optical  data  points  given  in  Table  2  are  shown,  normalized  to 
the  value  obtained  from  the  gas-metering  system,  in  Figs.  7  and  8.  Centerline  values 
obtained  from  both  the  gas-metering  system  and  from  probes  are  also  given  in  Table  2.  The 
spread  in  concentration  shown  in  Fig.  8  is  attributed  by  UTRC  to  uncertainty  in  the  mass 
spectrometer  measurements  (Ref.  4). 

The  temperature  profiles  for  the  FFB  were  measured  using  a  buttwelded  Ir/60%  lr/40% 
Rh  thermocouple  (0. 003-in. -dianr  wires)  coated  with  a  mixture  of  10%  beryllium  oxide  and 
90%  yttrium  oxide.  Radiation  corrections  were  made  using  the  thermal  balance  approach 
outlined  in  Ref.  4.  The  temperature  profiles  along  the  optical  path  for  the  FFB  are  shown  in 
Figs.  9  through  12.  The  radiation  correction  shown  on  Figs.  9  through  12  amounts  to  as 
much  as  200  K;  Ref.  4  estimates  the  error  in  the  correction  at  25  to  40  percent. 
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Tabie  2.  Operating  Conditions  for  Flat-Flame  Burner 


Data 

Point 

Number 

Temperature 

Profile 

Used 

Average 

Temperature, 

K 

Concentration 

Profile 

Used 

Centerline  NO 
Concentration,  ppm 
Metering  Probe 

Argon 

Mole 

Fraction 

Cal. 

HO 

Mole 

Fraction 

Cal. 

Distance 

From 

Burner 

Face, 

cm 

7 

16 

1,755 

Avg. 

1,390 

1,427 

0.5836 

0.414 

1.0 

10 

13 

913 

12 
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1 ,486 
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0.0678 

2.0 

13 

14 
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13 
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12 
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0.0680 

2.0 

17 

13 

913 

12 

686 

— 

0.8593 

0.0682 

2.0 

20 

14 

1,278 

10 

1,052 

1,119 
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— 
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Distance  from  Centerline,  cm 

Figure  7.  Normalized  concentration  profiles  of  the  flat-flame 
burner  at  room  temperature. 
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Figure  8.  Normalized  concentration  profiles  of  the  flat-flame 
burner  at  elevated  temperature. 
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Figure  9.  Temperature  profile  of  the  flat-flame  burner 
(1,000-K  setting). 


Figure  10.  Temperature  profile  of  the  flat-flame  burner 
(1,500-K  setting). 
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Distance  from  Centerline,  cm 

Figure  11.  Temperature  profile  of  the  flat-flame  burner 
(1,750-K  setting). 


3.2  AEDC  OPTICAL  SYSTEM  AND  MEASUREMENTS 

The  AEDC  UV  Resonance  Absorption  System  used  for  the  UTRC  measurements 
consisted  of  a  capillary  discharge  lamp  UV  source,  a  1/2-m  Jarrell-Ash  grating  spectrometer 
employing  a  Hamamatsu  R166UH  photomultiplier  detector,  associated  signal-conditioning 
electronics,  a  Hewlett-Packard  5480B  signal  averager,  and  a  Honeywell  5600C  fm  tape 
system;  the  system  is  shown  schematically  in  Fig.  13.  The  function  of  each  component  of  the 
system  is  described  in  more  detail  in  Refs.  1  through  3.  The  spectrometer  was  operated  at  a 
slit  width  such  that  the  half-width  of  the  triangular  slit  function  was  1.6  A.  Basic  to  the 
method  is  detailed  knowledge  of  the  fundamental  source  characteristics  —  line  positions  and 
intensities  and  the  temperature  in  the  discharge  region.  Reproducibility  of  the  source 
intensity  distribution,  as  shown  in  Fig.  14,  was  checked  at  UTRC  using  its  1.5-m  J-Y 
spectrometer  (Ref.  4).  UTRC’s  system  utilized  matched  optics  to  produce  a  parallel  source 
beam  1.0  cm  in  diameter  for  the  absorption  measurements. 
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Data  were  obtained  at  UTRC  by  scanning  the  spectrum  several  times  and  by  averaging 
the  spectrum  using  the  Hewlett-Packard  5480B  signal  averager.  The  averaged  spectrum  was 
then  recorded  on  the  fm  tape  and  on  an  X-Y  plotter.  A  reference  spectrum  was  obtained 
before  and  after  a  test  on  either  of  the  UTRC  calibration  sources,  and  several  data  points  at 
the  different  operating  conditions  given  in  Section  3.1  were  obtained.  Typical  spectral  data 
for  the  FGH  and  the  FFB  are  given  in  Figs.  15  and  16.  Data  for  the  absorption  cell*  were 
similar  and  are  not  shown  here. 


Problems  with  the  experimental  data  included  differences  in  reference  spectra  obtained 
before  and  after  a  test,  possibly  due  to  incomplete  purging  of  the  viewing  tube  and/or 
window  fogging.  Additional  differences  were  noted  between  hot  zero  (burner  on — no  seed 
of  NO)  and  cold  zero  (burner  off).  These  differences  contributed  to  the  optical  measurement 
uncertainty. 


Figure  12.  Temperature  profile  of 
the  flat-flame  burner 
(2,000-K  setting). 


*The  disagreement  of  these  data  with  the  original  AEDC  model  prompted  UTRC’s  investigation  of  the  model  and  the 
identification  of  errors  (Ref.  4).  The  subsequent  reexamination  of  the  model  by  AEDC  with  the  results  is  presented  herein. 
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Figure  13.  Schematic  of  AEDC  absorption  system  as  used  at  UTRC. 


Figurd  14.  Population  distribution 
of  excited  rotational 
states  of  the  A2  2  level 
of  NO  in  a  water-cooled 
discharge  tube  operated 
at  8  torr  with  2,800  v 
applied  and  containing 
12:3:1  mixture  (by 
volume)  of  A:N2:02. 


‘  45,000  46,000  47,000  48,000 

F(K'),  cm-1 


24 


A  EDC-TR-79-65 


■4 


Figure  15.  Typical  NO  (0,0)  7-band  absorption  spectra  from  flowing- 
gas  heater. 


Figure  16.  Typical  NO  (0,0)  7-band  absorption  spectra  from  flat- 
flame  burner. 
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3.3  DATA  TREATMENT 
3.3.1  Method  and  Procedure 

AEDC’s  original  plan  for  the  absorption  data  treatment  was  to  employ  the  iterative 
radial  inversion  model  of  the  UV  radiative  transfer  (described  in  Ref.  2)  for  a  single, 
homogeneous  zone.  In  that  model  the  NO  number  density,  or  concentration,  of  the 
homogeneous  zone  is  determined  by  successive  iterations  on  the  NO  density  until  the 
measured  and  calculated  second  bandhead  transmittances  are  the  same  to  within  some 
specified  small  increment.  Since  the  temperature  and  concentration  were  not  homogeneous 
along  the  FGH  or  FFB  optical  paths — as  shown  by  the  measured  profiles  given  in  section 
3.1— this  plan  was  dropped.  Instead,  it  appeared  that  the  best  approach  to  a  high- 
temperature  calibration  of  the  model  was  to  compare  transmissivities  calculated  using  the 
measured  temperature  and  NO  concentration  profiles  and  to  break  the  path  up  into  several 
homogeneous  zones  for  which  averaged  temperatures  and  NO  concentrations  had  been 
assigned. 

The  radiative  transfer  of  line  radiation  along  nonhomogeneous  paths  is  complicated  by 
the  changing  line  shape  along  the  path.  This  problem  is  treated  in  several  classic  references, 
including  Goody  (Ref.  16).  Beer’s  law,  which  is  sometimes  used  erroneously  for  the  radiative 
transfer,  state  that  the  monochromatic  absorption  coefficients  of  zones  of  different 
properties  are  additive.  Sometimes  the  integrated  line  transmittances  are  treated, 
erroneously,  as  the  product  of  transmittances  of  the  different  zones.  This  procedure  is  valid 
only  if  the  source  line  is  much  narrower  than  the  absorption  lines. 

The  procedure  used  here  for  the  radiative  transfer  calculation  is  illustrated  in  Fig.  17  and 
is  detailed  for  a  single  line,  j,  by 


L 

l 

t  .(id  =  I°0  (i /)  exp 

£~aX(! 

~r"a 

r^: 

•b 

r^j 

! 

3  j 

r=i 

(15) 


TiWd^ 


(16) 


Expressed  in  words,  Eqs.  (15)  and  (16)  state  that  the  intensity  of  the  transmitted  line,  tj,  at 
the  right  (receiver)  side  of  the  nonhomogeneous  path  made  up  of  zones  of  constant 
properties  is  calculated  by  starting  the  calculation  on  the  left  (source)  side  (see  Fig.  17)  and 
determining  the  frequency  function  in  Eq.  15  for  the  first  zone  by  breaking  the  line  shape  up 
into  a  sufficient  number  of  monochromatic  points  to  define  the  line  shape  (e.g.,  100 
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Figure  17.  Illustration  of  spectral  line  radiative  transfer  through 
nongomogeneous  media. 


frequency  points  were  used  for  the  calculations  reported).  The  result  is  a  new  line  profile, 
to  be  used  for  zone  2,  in  which  the  absorption  coefficients  for  all  the  lines,  i, 
contributing  to  absorption  of  source  line,  j,  are  determined  by  the  temperature,  T2, 
pressure,  p2,  and  NO  density,  n2,  in  zone  2.  The  result  of  the  zone  2  calculation  is  the  source 
line  profile  to  be  used  for  zone  3;  the  zone  3  calculation  is  used  for  zone  4  and  so  forth.  At 
the  right  boundary  of  the  path,  the  final  shape  of  the  transmitted  profile,  tj(j'),  has  been 
determined,  and  the  transmitted  intensity  integral,  t  j  [Eq.  (16)],  is  evaluated.  The  integral 
transmission  values  are  stored  for  each  line;  the  computer  program  that  sums  the 
contributions  of  each  line  within  the  instrument  bandpass  is  then  run.  The  line  shapes  shown 
in  Fig.  17  are  exaggerated  to  emphasize  the  problem.  The  shapes  represent  the  effect  of  a 
cold  boundary  with  a  line  width  smaller  than  the  source  lien  width  and  with  subsequent 
broad  lines  corresponding  to  a  high-temperature  path.  Such  a  situation  would  be 
encountered  for  a  jet  engine  exhaust  under  altitude  conditions. 


The  revised  AEDC  model  described  in  Section  2.0  was  used  along  with  the  above- 
described  radiative  transfer  model  to  calculate  transmissivities  at  the  second  bandhead 
expected  for  the  FGH  and  FFB  property  profiles  given  in  Section  3.1.  The  23  zones  used  are 
shown  in  Fig.  18.  Local  values  of  the  broadening  parameter,  (J,  were  computed  using  Eq. 
(14)  to  determine  the  broadening  coefficient. 
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Figure  18.  Diagram  of  zone  widths  for  computer  calculation 
of  transmission  for  the  flowing-gas  heater  and 
flat-flame  burner. 


3.3.2  Results 

The  computed  values  of  the  second  bandhead  transmittance  for  the  absorption  cell,  the 
FGH,  and  the  FFB  are  given  in  Tables  3,  4,  and  5,  respectively;  the  procedure  described  in 
the  preceding  paragraph  was  used  with  the  temperature  and  NO  concentration  profiles  given 
in  Section  3.1.  The  centerline  values  of  concentration  listed  in  Tables  1  and  2  were  used,  and 
a  uniform  pressure  of  1  atm  was  assumed.  For  several  of  the  data  points,  profiles  as 
measured  by  probes  were  not  supplied,  and  the  profiles  shown  by  the  solid  curves  in  Figs.  4 
through  12  were  used. 

The  transformation  of  the  differences  in  measured  and  calculated  second  bandhead 
transmissivity  to  expected  differences  in  NO  concentration  was  computed  by  the  procedure 
suggested  by  UTRC  (Ref.  4),  using  the  following  equation: 


percent 


[NO]  difference 


(  in  t 

[  _ meat 

\  £n  r  . 

\  cal 


100 


(17) 


The  results  of  this  calculation  are  also  shown  in  Tables  3,  4,  and  5.  Such  an  expression  may 
serve  as  a  basis  for  comparison  but  cannot  be  interpreted  as  a  very  good  estimate  of  error 
since  it  applies  only  to  single  lines  in  the  narrow  source  line  limit  and  to  a  homogeneous 
path.  Another  approach  was  taken  by  AEDC  for  a  few  data  points,  the  results  of  which  are 
shown  in  Table  5.  In  this  approach,  the  centerline  concentration  was  varied  and  the 
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transmittance  calculation  performed  several  times  until  a  fit  with  the  measured 
transmittance  was  effected.  This  AEDC  method  gave  differences  in  predicted  and  measured 
transmittance  usually  about  five  percent  less  than  those  obtained  using  Eq.  (17). 


Table  3.  Results  of  Optical  Transmission  Measurements  Made  Through 
Absorption  Cell  (293  K) 


Data 

Point 

Number 

Pressure, 

atm 

NO 

Concentration, 

c,  Ppm 

Transmissivity, 
t,  measured 

Transmissivity , 
t,  model 

At* 

** 

Ac  , 
percent 

i 

1 .0 

458 

0.482 

0.454 

0.028 

“7.5 

2 

1.0 

93 

0.835 

0.842 

0.007 

4;  8 

3 

0.0395 

2,080 

0.673 

0.627 

0.046 

-15.1 

4 

1.0 

250 

0.659 

0,637 

0.022 

-7.5 

*Al  =  T  ,  -  -  T 

model  meas 

(£n  t  \ 

meas  l 

- 1  I  100  percent 

fn  Tmodel  / 


Table  4.  Results  of  Optical  Transmission  Measurements  Made  Through 
Flowing-Gas  Heater 


Data 

Point 

Number 

Centerline  NO 
Concentration, 
c ,  ppm 

Transmissivi ty , 
t,  measured 

Transmissivity , 
t,  model 

ft 

At 

ftft 

Ac  , 
percent 

2 

115 

0.836 

0.826 

-0.010 

-6.3 

8 

115 

0.862 

0.834 

-0.028 

-18.2 

8 

112  (p) 

0.862 

0.838 

-0.024 

-16.0 

9 

301 

0.686 

0.634 

-0.052 

-17.3 

9 

287  (p) 

0.686 

0.646 

-0.640 

-13.8 

11 

115 

0.881 

0.906 

0.025 

28.4 

12 

301 

0.743 

0.775 

0.032 

16.5 

14 

1 15 

0.898 

0.908 

0.010 

11.5 

15 

301 

0.  748 

0.  776 

0.028 

14.5 

17 

301 

0.  731 

0.  762 

0.031 

15.3 

19 

115 

0,952 

0.953 

0.001 

2.2 

19 

105  (p) 

0.952 

0.957 

0.003 

11.9 

20 

301 

0.  879 

0.876 

-0.003 

-2.6 

20 

280  (p) 

0.879 

0.884 

0.005 

3.7 

23 

115 

0.923 

0.947 

0.024 

47.1 

24 

301 

0.840 

0.869 

0.029 

24.2 

27 

535 

0.762 

0.793 

0.031 

17.2 

27 

520  (p) 

0.  762 

0.  798 

0.036 

20.5 

28 

448 

0.812 

0.831 

0.019 

12.5 

28 

430  (p) 

0.812 

0.837 

0.025 

14.7 

(p)  -  probe  value 


Tmodel  Tmeas 
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Table  5.  Results  of  Optical  Transmission  Measurements  Made  Through 
Flat- Flame  Burner. 


Data 

Centerline  NO 

Point 

Concentration, 

Transmissivity 

Number 

c,  ppm 

t,  measured 

7 

1 ,390 

0.824 

7 

1,427  (p) 

0.824 

10 

1,394 

0.607 

10 

1,486  (p) 

0.607 

13 

1,362 

0.728 

13 

1,459 

0.728 

16 

1 ,076 

0.674 

17 

686 

0.  790 

20 

1,052 

0.  746 

20 

1,119  (p) 

0.746 

21 

670 

0.832 

24 

1,000 

0.  817 

24 

1,107  (p) 

0.817 

25 

1,141 

0.792 

26 

750 

0.861 

27 

607 

0.896 

30 

1 ,584 

0.779 

30 

1.671  (p) 

0.779 

33 

301 

0.578 

34 

214 

0.638 

35 

115 

0.798 

(p)  -  probe  value 


■k  + 

Transmissivity,  A  Ac  ,  Ac^ 


model 

At 

percent 

percent 

0.854 

0.030 

22.7 

0.851 

0.027 

20.0 

— 

0.587 

-0.020 

-6.3 

— 

0.567 

-0.040 

-12.0 

— 

0.  759 

0.031 

14.6 

— 

0.744 

0.016 

7.4 

— 

0.661 

-0.014 

-4.7 

— 

0.  766 

-0.024 

-11.6 

— 

0.807 

0.061 

36.7 

— 

0.  796 

0.050 

28.4 

— 

0.872 

0.040 

34.3 

— 

0.851 

0.034 

25.4 

23.0 

0.837 

0.020 

13.4 

12.9 

0.832 

0.040 

26.9 

23.5 

0.886 

0.025 

23.4 

20.6 

0.906 

0.010 

11.2 

7.6 

0.836 

0.057 

39.4 

— 

0.828 

0.049 

32.  3 

— 

0.583 

0.005 

1.6 

— 

0.676 

0.038 

14.  8 

— 

0.806 

0.008 

4.0 

— 

5fAi  =  t  ,  1  —  t 

model  meas 

An  t  \ 

/  meas  \ 

Ac  =1 - 1  1  100  percent 

' Tmodel  / 


t  determined  by  matching  centerline  concentration  to  fit  model. 


3.3.3  Discussion 


The  uncertainties  introduced  by  the  nonhomogeneous  profiles  of  temperature  and 
concentration  in  the  FGH  and  FFB  into  the  “calibration”  of  the  NO  UV  absorption  method 
at  high  temperatures  made  the  calibration  somewhat  less  than  ideal.  In  fact,  by  judicious 
choice  of  profiles  on  most  of  the  data  points,  the  measured  and  calculated  transmittances 
could  be  made  to  agree.  For  example,  the  uncertainty  in  the  radiation  correction  to  the 
temperature  for  the  FFB  always  permits  the  transmissivity  to  be  matched  by  suitable  choice 
of  temperature  within  the  error  bounds.  It  should  be  noted  that  the  preponderence  of  the 
data  indicates  an  over-prediction  of  the  NO  concentration  (over-prediction  of  the 
transmittance).  Taken  overall,  the  root-mean-square  average  of  the  departure  of  the 
calculated  NO  concentration  using  Eq.  (16)  is  as  follows: 

Absorption  cell  9  percent 
Flowing-gas  heater  18  percent 
Flat-flame  burner  21  percent 

These  results  are  considered  satisfactory,  considering  the  uncertainties  in  the  profiles  of 
temperature  and  concentration. 
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Although  it  represents  the  correct  physics,  the  radiative  transfer  calculational  procedure 
introduced  here  would  not  be  expected  to  be  very  different  from  a  Beer’s  law  treatment  since 
the  absorption  line  widths  in  these  experiments  were  considerably  larger  than  the  source  line 
width  (i.e.,  from  about  2  to  4).  However,  for  applications  to  turbine  engine  exhausts  at 
altitude  conditions,  the  Beer’s  law  procedure  would  be  expected  to  be  considerably  in  error. 


4.0  SUMMARY  OF  RESULTS 

The  results  obtained  by  AEDC  in  Phase  1  of  the  Interagency  NO  Measurement 
Investigation  are  summarized  as  follows: 

1.  Transmissivity  measurements  of  the  NO  (0,0)  y  band  were  made  at  the  United 
Technology  Research  Center  laboratories  using  the  calibration  cells — room- 
temperature  absorption  cell,  flowing-gas  heater  (293  to  833  K),  and  flat-flame 
burner  (230  to  1,755  K).  NO  concentration  and  temperature  profile  data 
obtained  on  the  flowing-gas  heater  and  flat-flame  burner  by  UTRC  using  probe 
methods  were  supplied.  These  data  indicated  nonhomogeneities  which 
complicated  the  calibration  procedure  established  for  a  homogeneous  path. 

2.  UTRC  personnel  discovered  errors  in  the  NO  resonance  line,  radiative  transfer 
model  developed  at  AEDC;  AEDC  made  the  necessary  revisions,  which  required 
a  redetermination  of  parameters.  The  completed  model  redevelopment  produced 
a  much-improved  model  and  more  reliable  parameters.  (Independently,  UTRC 
also  developed  a  revised  model,  different  in  detail  but  substantially  the  same  as 
the  revised  AEDC  model,  and  also  acquired  moderately  high  resolution  line  data 
and  determined  an  independent  set  of  parameters.) 

3.  The  newly  obtained  broadening  cross  sections  are  approximately  equal  for  N2, 

A,  and  C02  (about  547  ±  53  x  10  16  cm2).  Broadening  by  H20  (for  which 
experimental  data  are  difficult  to  obtain)  and  02  were  assumed  to  have  the  same 
cross  section  as  N2,  A,  and  C02.  This  is  important  for  the  H2/02  flat-flame 
burner  used  by  UTRC  in  which  as  much  as  40  percent  of  the  effluent  volume  is 
H20.  This  assumption  results  in  added  uncertainty  in  the  calculated 
transmissivities. 

4.  A  radiative  transfer  calculation  procedure  for  spectral  lines  transmitted  through 
nonhomogeneous  media  was  developed  and  applied  to  the  data  obtained  at 
UTRC  by  AEDC.  The  procedure  avoids  incorrect  applications  in  the 
assumptions  of  Beer’s  law,  which  are  reasonable  for  this  data  where  the  absorber 
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line  half-widths  are  considerably  greater  than  the  source  line  half-widths.  Serious 
difficulty  with  the  Beer’s  law  assumption  would  be  expected  for  sources  such  as 
turbine  engine  exhausts  at  simulated  altitudes. 

5.  The  calibration  results  for  the  AEDC  system  from  Phase  1  of  the  Interagency 
program,  though  not  ideal  because  of  the  nonhomogeneity  of  the  paths  and 
because  of  the  uncertain  temperature  and  concentration  profiles  in  the 
calibration  devices,  were  deemed  satisfactory.  An  average  prediction  of  about 
20-percent  uncertainty  in  NO  concentration  can  be  made  using  the  AEDC  UV 
resonance  line  absorption  technique,  based  upon  these  results. 

REFERENCES 

1.  Few,  J.  D.,  Bryson,  R.  .1.,  and  McGregor,  W.  K.  “Evaluation  of  Probe  Sampling  versus 
Optical  In  Situ  Measurements  of  Nitric  Oxide  Concentrations  in  a  Jet  Engine 
Combustor  Exhaust.”  AEDC-TR-76-180  (ADA034726),  September  1971. 

2.  Few,  J.  D.,  McGregor,  W.  K.,  and  Glassman,  H.  N.  “Comparison  of  UV  Absorption 
Measurements  with  Probe-Sampling  Measurements  of  Nitric  Oxide  Concentrations  in  a 
Jet  Engine  Combustor  Exhaust.”  AEDC-TR-76-134  (ADA030005),  September  1976. 

3.  Few,  J.  D.  “Optical  Measurement  of  NO  and  NO2  in  the  Exhaust  of  an  F101-GE-100 
Engine  at  Simulated  Altitudes.”  AEDC-TR-77-75  (AD-A047882),  December  1977. 

4.  Dodge,  L.  G,,  Colket,  M.  B.,  Ill,  Zabielski,  M.  F.,  Dusek,  J.,  and  Seery,  D.  J.  “Nitric 
Oxide  Measurement  Study:  Optical  Calibration.”  DOT-FAA  Contract  No. 
FA77WA-4081 .  (In  review.) 

5.  Davis,  M.  G.,  McGregor,  W.  K.,  Few,  J.  D.,  and  Glassman,  H.  N.  “Transmission  of 
Doppler  Broadened  Resonance  Radiation  Through  Absorbing  Media  with  Combined 
Doppler  and  Pressure  Broadening  (Nitric  Oxide  7-Bands  as  an  Example).”  AEDC- 
TR-76-12  (ADA021061),  February  1976. 

6.  Herzberg,  G.  Spectra  of  Diatomic  Molecules.  D.  Van  Nostrand  Co.,  New  York,  1950. 

7.  Tatum,  J.  B.  Astrophysical  Journal.  Supplement  Series.  Vol.  14,  Supplement  No.  124, 
1967. 

8.  Pery-Thorne,  A.  and  Banfield,  F.  G.  “Absolute  Oscillator  Strength  of  the  (0,0)  Band 
of  the  Gamma  System  of  Nitric  Oxide  by  the  Hook  Method.”  Journal  of  Physics.  B: 
Atomic  and  Molecular  Physics,  Vol.  3,  No.  7,  July  1970,  pp.  1011-1019. 


32 


A  EDC-TR-79-65 


9.  Keefer,  D.  R.  “Level  Populations  for  Molecules  Intermediate  to  Hund’s  Cases  (a)  and 
(b):  The  NO  X2II  State  as  an  Example.”  (Submitted  to  Journal  of  Quantitative 
Spectroscopy  and  Radiative  Transfer.) 

10.  Penner,  S.  S.  Quantitative  Molecular  Spectroscopy  and  Gas  Emissivities.  Addison- 
Wesley  Co.,  Reading,  Mass.,  1959. 

11.  Deeszi,  I.  Acta  Physica.  Vol.  9 

12.  Engleman,  R.,  Jr.,  Rouse,  P.  E.,  Peck,  H.  M.,  and  Baiamonte,  V.  D.  “Beta  and 
Gamma  Systems  for  Nitric  Oxide.”  Los  Alamos  Report  LA-4364,  Los  Alamos 
Scientific  Laboratory,  University  of  California,  New  Mexico,  July  1970. 

13.  Hasson,  V.,  Farmer,  A.  .1.  D.,  Nicholls,  R.  W.,  and  Anketell,  J.  “Application  of 
Dispersion  Techniques  to  Molecular  Band  Intensity  Measurements.  II:  Oscillator 
Strength  of  the  (0,0)  Band  of  NO-y  (A2£-X3n)  System.”  Journal  of  Physics.  B: 
Atomic  and  Molecular  Physics,  Vol.  5,  No.  6,  June  1972,  pp.  1248-1254. 

14.  Farmer,  A.  J.  D.,  Hasson,  V.,  and  Nicholls,  R.  W.  “Absolute  Oscillator  Strength 
Measurements  of  the  {v"  =  0,  v'  =  0-3)  Bands  of  the  (A2£-X2n)  y-System  of  Nitric 
Oxide.”  Journal  of  Quantitative  Spectroscopy  and  Radiative  Transfer,  Vol.  12,  No.  4, 
April  1972,  pp.  627-633. 

15.  Benson,  R.  and  Samuelson,  G.  S.  “Oxides  of  Nitrogen  Transformation  While  Sampling 
Combustion  Products  Containing  Carbon  Monoxide  and  Hydrogen.”  Presented  at 
Combustion  Institute,  Fall  Meeting,  LaJolIa,  California,  October  18-20,  1976.  Paper 
76-39. 

16.  Goody,  R.  M.  Atmospheric  Radiation.  I:  Theoretical  Basis.  Clarendon  Press,  Oxford, 
1964. 

17.  Mitchell,  A.  C.  G.  and  Zemansky,  M.  W.  Resonance  Radiation  and  Excited  Atoms. 
Cambridge  University  Press,  London,  1971. 

18.  Present,  R.  D.  Kinetic  Theory  of  Gases.  McGraw-Hill,  New  York,  1958. 


33 


A  EDC-TR  -79-65 


APPENDIX  A 

DERIVATION  OE  COLLISION AL  BROADENING  RELATIONS 

Some  confusion  exists  in  the  literature  in  regard  to  the  origin  of  expressions  for 
collisional  “half-widths,”  “optical  collision  diameters,”  and  “optical  cross  sections.”  This 
appendix  is  intended  to  clarify  these  expressions  as  they  are  used  in  this  report. 

The  full  width  at  half-maximum,  (Ai>)l,  of  a  collision-broadened  spectral  line  is  defined 
as  Zl/7t,  where  ZL  is  the  number  of  collisions  per  unit  time  per  atom  giving  rise  to  the 
broadening  of  a  line  (Ref.  17).  The  Lorentz  frequency  function  is  classically  used  to  describe 
the  shape,  bL(v),  of  a  natural  or  collision-broadened  line, 


bL(v) 


constant  {Av)^2 

(A^)l  2  +  (i/  —  )2 


(A-l) 


where  v°  is  the  frequency  at  line  center.  Now,  if  the  concepts  of  kinetic  theory  are  used,  the 
mean  collision  frequency  between  two  molecules,  a  and  b,  per  number  of  molecules,  a,  is 
given  by 


[(/V)l]  ^  =  Zab  =  (77dab)nbvab 


(A-2) 


where  nb  is  the  number  density  of  molecules,  b;  vab  is  the  mean  speed  between  the  molecules; 
and  dab  is  the  diameter  of  the  sphere  of  influence  between  the  molecules  involved  in  the 
collision.  Thus,  7rd lb  is  an  area  which  may  be  designated  as  a  “cross  section”  for  the 
collision,  aab.  The  “optical  collisional  diameter”  for  broadening  is  defined  as  dab.  Both 
terms  are  in  reality  unknown  and  are  directly  related  to  the  probability  that  a  collision 
promotes  a  departure  from  the  stationary  state  probability  of  transition.  The  mean  speed  for 
molecules  having  a  Maxwell-Boltzmann  velocity  distribution  function  is  (Ref,  18) 


v 


ab 


(A-3) 
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where  x  is  Boltzmann’s  constant,  T  is  the  temperature,  Mab  is  the  the  reduced  mass,  and  Ma 
and  Mb  are  the  masses  of  the  molecules  a  and  b,  respectively.  The  density  of  molecules,  b,  is 
given  by 


n 


b 


(A-4) 


where  p  is  the  pressure  of  the  gas  mixture  and  fb  is  the  mole  fraction  of  molecules,  b. 
Combining  Eqs.  (A-2)  through  (A-4)  results  in  the  following: 


(A-5) 


If  more  than  one  molecule  b  contributes  to  the  broadening  of  spectral  lines  emitted  or 
absorbed  by  molecule  a,  then 


(A-6) 


If  the  lines  are  also  broadened  by  the  Doppler  effect,  the  line  shape  is  given  by  the  Voigt 
profile  (Ref.  17): 


bvd 


i 

77 


ffe-r2 

Q2  +  (co-  y)2 


dy 


(A-7) 


where  y  is  a  variable  of  integration,  w  is  the  Doppler  function 


Cl) 


2  Vfn  2 


(ACt 


(A-8) 
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with  v°  the  line  center  frequency  and  (ArOD  is  the  full  width  at  half-maximum  of  the  Doppler 
line  shape  function.  The  latter  is  given  by 


bDM  =  e-"2 


(A-9) 


and  the  Doppler  half-width  is  (Ref.  17) 


2  fn  2  kT 
M  c2 


(A- 10) 


where  c  is  the  velocity  of  light.  The  parameter  (J  in  Eq.  (A-7)  is  the  “collisional  broadening 
parameter”  and  is  defined  as 


a 


(At 


(At>)r 


(A-l  1) 


Now,  inserting  Eqs.  (A-6)  and  (A-10)  into  (A- 11)  and  putting  v°  =  c/X  where  X  is  the 
wavelength  of  the  line  give  the  following: 


a 


.3./  2 


(A- 12) 


or 


a  = 


(A- 13) 
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where  C  is  the  “collision-broadening  coefficient”  defined  as 


c  =  -171  t  th<Jabl1  +  K 

KTT  b  b 


(A- 14) 


In  terms  of  the  “optical  collision  diameter”  Eq.  (A-14)  becomes 


!-\,n  b 


(A- 15) 
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